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PROPAGATION  OP  ELECTROMAGNETIC  WAVES  IN  PARAMETRIC  MEDIA  IN  THE 
APPROXIMATION  OF  GEOMETRIC  OPTICS. 

V.  N.  Krasil'nikov. 

§1.  Posing  of  the  problen. 

During  the  passage  of  waves  through  media  whose  parameters  are 
changed  in  time  (subseguently  let  us  call  them  simply  parametric), 
occur  special  effects,  for  example  a change  in  the  spectral 
composition  of  signal.  Therefore  to  their  study  traditional  approach 
[1]  proves  to  be  inapplicable,  and  is  required  special  examination, 
in  order  to  explain  the  possibility  of  using  the  laws  of  geometric 
optics.  In  this  case,  systematically  it  is  more  justified  to  analyze 
directly  the  equations  of  Maxwell,  but  not  to  reduce  them  for 
different  special  cases  to  one  differential  equation  for  cortain 
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scalar  function,  as  this  frequently  is  done  [1]. 

Since  the  approxiaat ion  of  geoaetric  optics  knowingly  does  not 
befit  for  the  description  of  fields  in  the  vicinity  of  their  sources, 
the  usual  fornulation  of  the  problen  of  geoaetric  approximation 
consists  of  finding  of  wave  process  of  space,  if  is  known  the  value 
of  field  on  certain  initial  surface  Sq.  The  latter  can  be  selected  in 
the  wave  zone  of  eaitter.  It  can  pass  in  those  regions  of  space, 
where  there  is  still  no  perturbation  action  of  heterogeneous  and 
paraaetric  properties  of  the  aediua.  Surface  Sq  is  convenient  +o 
coabine  with  the  position  of  wave  front  (in  aonochromatic  field  ^his 
is  the  constant'phase  surface)  . 

Outside  surface  Sq  let  us  consider  the  aedium  isotropic, 
heterogeneous  and  paraaetric;  its  properties  are  describe!  by  three 
functions:  dielectric  and  aagnetic  constants  •(•»,  r,  t)  and  p(«,  r,  t)‘  and 
conductivity  •»)(»,  r,  f)-  If  priaary  radiation  (field  on  surface  Sq)  can 
be  considered  aonochroaatic  (frequency  wq) , then  due  to  the 
paraaetric  nature  of  the  aediua  the  spectral  coaposition  of  sional 
will  be  changed  in  certain  interval  of  frequencies  Au  near 


Page  77. 


As  we  shall  see  further  on,  relative  frequency  drift  in  all  cases  of 
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practical  interest  is  saall,  and  it  is  reasonable  to  disregard  the 
dispersion  of  the  characteristics  of  the  sediuo  in  interval  Au , i.e., 
to  take  the  values  of  functions  and  when  u = uo«  Then  in  the 

Gaussian  systen  of  units  the  equations  of  Naxwell  it  is  possible  to 


write  thus: 


rotE+-J-4-(pH)  = 0. 

rotH-4-4(.E)  = ip-E. 


(») 


Let  us  note  that  on  the  strength  of  the  linearity  of  system  ( 1) 
and  for  parametric  media  it  is  possible  to  utilize  principle  of 
superposition  and  to  describe  the  dependence  of  electromagnetic 
vectors  E-  and  M (but  not  the  parameters  ^ ,y~  and  !)  on  time 

with  the  aid  of  thecomposite  factor  exp{-iuQt}  - we  will  consider 
that  on  the  initial  surface  Sq  ve  have  precisely  this  time  behavior 
of  field. 


! 


Let  us  impose  on  the  properties  of  the  medium  limitations  which 
make  it  possible  to  apply  to  the  analysis  of  the  process  of  the 
propagation  of  waves  the  method  of  geometric  optics. 


Pirst,  it  is 
heterogeneous,  in 
be  changed  little 
(at  a given  place, 
values: 


well  known  [1]  that  the  medium  must  be  weakly 
other  words,  its  characteristics  £,p.and  fj  must 
at  distances  of  the  order  of  the  local  wavelength  x 
at  a given  instant),  i.e.,  if  L is  smallest  of  thv 


/ l«l  [ li«l  / _ 1 1 1 
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then 

-r«>-  (2) 

Further  let  us  consider  the  medium  weakly-parametric,  i.  e.  , assume 
that  its  properties  little  are  changed  during  a time  on  the  order  of 
M-»,  where  u is  local  frequency  of  field.  This  condition  takes  the 
form 

(3) 

where  T is  smallest  of  the  times: 


7'.  = | 

1 ( 

dt 

I* 

. 7-,= 

1 

~Si 

•ST 

dt 

Furthermore,  it  is  reasonable  to  consider  that  the  attenuation 
of  wawes  at  a distance  of  the  order  of  wavelength  is  slight 
(otherwise  the  very  problem  of  propagation  does  not  appear).  Then  it 
is  possible  to  introduce  the  third  low  parameter  of  the  problem 

^«1.  (4) 

characterizing  the  relative  role  of  the  bias  currents  and  conduction 
currents  (latter  are  small).  Condition  (4)  is  not  in  principle 
necessary,  but  it  makes  it  possible  to  obtain  the  results  in  a more 
compact  and  visible  form. 


Sometimes  inequalities  (3)  and  (4)  are  conveniently  written  in 
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the  fora  of  the  ratio  of  wavelength  X to  certain  large  graphic  scale. 
So,  if  V - the  velocity  of  propagation  of  electromagnetic  wave,  then 
(3)  is  equivalent  to  the  inequality 


where  Lt — while  (4)  it  is  possible  to  present  thus: 

-^«l.  (4a> 

'■’I 

determines  the  distance  at  which  the  wave  amplitude  decreases 
due  to  losses  in  the  medium  by  a factor  of  e. 


Let  us  search  for  the  solution  of  system  of  equations  (1)  in  the 


form 


E(r,  0 = Eo('’.  0 (5) 

H(r.  /)  = H,(r, 

where  are  clearly  isolated  the  common  for  fields  ^ and  ^ rapid ly 
changing  phase  factor  and  the  slowly  changing  amplitude  vectors 

Ht  • assume  that  from  the  existence  of  the  small  parameter: 

(2),  (3)  and  (4)  ensues  the  possibility  of  the  separation  of  all 

entering  the  equations  expressions  into  two  categories: 


^ 1)  the  values  of  basic  order;  to  them  are  related 

■Jf  ■’ 


2)  the  small  first-order  quantity,  at  least  one  of  out  small 
parameters.  This:  rotE,.  rot  H,;  vE,:  E,; 
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At  the  basis  of  this  classification  lie  evaluations  of  the  type 

and  obvious  deterainations  for  local  frequency  “=1^  local 

wavelength  X = r?l  Let  us  emphasize  that  we  still  must  show  that  the 

iv?r 

taken  separation  of  values  into  orders  actually  follows  from  the 
restrictions  placed  on  the  properties  of  the  medium. 


§2.  Phase  relationships 


Let  us  apply  the  method  of  successive  approximations  to  t h*^- 
system  of  equations  of  Maxwell. 

Page  79, 


Consideration  of  the  terms  of  basic  order  leads  to  the  relationships 

X H,  = 7 

whence  for  phase  function  « follows  the  generalized  equation  of  the 


(6) 


eikonal 
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physical  sense  of  which  is  very  simple;  it  describes  the  propagation 
of  constant-phase  surface  * = const  at  the  local  speed  = K ep  in 
directions  ±7?- it  is  not  difficult  to  interpret  relationships  (6): 
vectors  and  are  oriented  orthogonally  to  each  other  and  to 

the  direction  of  propagation  of  wave 

FOOTNOTE  1.  Here  and  subsequently  we  for  definiteness  will  analyze 
the  waves  which  go  in  direction  tvvp  .endfootnote. 


Between  their  values  there  is  connection  through  the  local  impedance 

1=:-SL 

of  the  wediuB.  Introducing  unit  vector  in  direction  of  beam 
from  equations  (6)  and  (7)  we  have 

Eo=/fllXH„).  (8) 

In  other  words,  with  an  accuracy  to  the  small  first-order 
quantities  the  local  structure  of  the  field  of  electromagnetic  wave 
coincides  with  the  same  for  the  case  of  plane  wave  in  the  medium  with 
constants  £ and 

We  already  mentioned  that  typical  for  geometric-optical 
calculations  is  this  posing  of  the  question,  when  field  is  considered 
assigned  at  certain  phase  wave  front  Sq- 


For  the  slightly 
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inhomogeneous  and  weak ly- parametric  media  the  field  of  concentrated 
source  at  distances  of  several  wavelengths  from  it  will  still  li*tle 
differ  from  the  field  which  it  would  create  in  the  homogeneous  medium 
with  the  time-constant  parameters.  Therefore  with  sufficient  for 
practical  calculations  accuracy,  the  field  in  the  vicinity  of  emitter 
can  be  considered  known.  In  this  case,  it  is  obvious  that  in  the  near 
7one  of  the  source  where  play  role  static  and  induction  fields, 
condition  (8)  does  not  obtain;  this  corresponds  to  the  known  fact  of 
the  inapplicability  of  the  geometric-optical  approximation  for  the 
calculation  of  the  local  fields  of  emitter. 

Thus,  let  the  constant-phase  surface  lie  in  the  wave  zone  of 
source,  and  at  certain  moment  of  time  tp  the  value  of  the  phase  in  Sq 
is  #0  ~ ♦(fo)-  I'tom  the  generalized  equation  of  the  eikonal  (7)  it 
follows  that  this  phase  front  propagates  at  a rate  of  v in  the 
direction  of  the  normal  to  itself  (the  latter,  by  definition,  is  the 
direction  of  beam)  . 

Page  80. 

Having  in  principle  radiation  picture,  we  can  trace  in  detail  ♦he 
motion  of  this  phase  of  wave  # = ♦(,  in  time.  In  this  case,  the 
inconstancy  in  time  of  speed  forces  us  during  calculation  of 

the  phase  picture,  which  corresponds  to  Pig.  1,  to  take  at  each  point 
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of  space  the  valua  of  velocity  at  that  noment  of  time  when  the  phase 
front  in  question  passes  through  this  point.  Let  us  agree  this 
delaying  value  of  the  velocity  subsequently  to  designate  [v].  It  is 
obvious  that  the  spatial  distribution  of  [v]  in  parametric  media  is  a 
function  of  the  starting  time  to,  and,  generally  speaking,  the 
radiation  picture  of  the  transfer  of  another  phase  of  the  wave  q, 
"starting"  from  the  initial  surface  So  at  another  moment  of  time  t'o, 
will  be  different.  Onder  such  unsteady  conditions  it  is  convenient  to 
construct  distributions  of  the  beam,  recording  not  current  time  t, 
but  the  starting  time  to  in  the  manner  shown  in  Fig.  1. 


It  is  easy  to  see  that  propagation  time  of  this  phase  of  wave 
from  point  Mq  in  initial  surface  to  certain  point  M,  which  lies  on 
the  beam,  which  left  at  the  moment  of  time  tn.  is  ecjual 


M 


where  the  integral  is  taken  on  the 
phase  front,  i.e.,  the  "length"  of 


I 


where  the  function  [v]  also  must  b 
of  beam. 


(9) 

beam,  and  the  path,  passed  by 
beam, 

(10) 

taken  at  the  appropriate  points 


During  the  calculation  of  beam  trajectories,  can  prove  to  be 
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useful  the  differential  equation  of  beam.  For  its  derivation  let  us 
examine  at  moment  of  time  t so  small  a section  of  phase  front  « = 
that  it  would  be  possible  to  consider  it  plane.  After  differen^-ially 
short  time  dt,  phase  front  will  be  shifted  in  space  in  direction 
in  distance  ds  = vdt  (Fig.  2). 


DOC  = 77218400 


PAGE 


ft 


Pig.  1. 

Page  81. 

If  speed  V (t)  is  not  identical  for  all  points  of  front,  then  its 
position  at  moment  of  time  t + dt  will  not  be  parallel  to  the  initial 
one.  Being  limited  to  the  small  first-order  quantities,  wa  will 
obtain  the  rotation  of  phase  front  (and,  this  means,  the  unit  vector 
of  beam  I ^ through  certain  angle  da: 

V = l + dl,  ldH  = da, 


where  rflil  is  caused  by  difference  from  zero  of  gradient  of  speed  v in 
the  direction  orthogonal  to  Simplest  geometric  calculation  gives 
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where  = — is  a transverse  to  1 part  ^v.  Introducing 
refractive  index  n = c/v,  we  come  to  the  following  differential 
equation  of  bean: 

/,(«!)=  (11) 

It  is  interesting  to  note  that  precisely  to  the  same  equation 
are  subordinated  the  beams  in  the  medium  with  constant  parameters 

[2] ,  but  in  our  case,  as  can  be  seen  from  the  given  derivation,  it  i 
correct  only  during  the  analysis  of  the  family  (or  distribution)  of 
beams  united  by  the  total  starting  tine  to. 

In  concluding  this  paragraph,  recall  that  sometimes  during  the 
construction  of  the  equations  of  geometric  optics  as  the  basis  of 
analysis  is  used  the  postulate  about  the  extreme  character  of  beams 

[3]  - the  so-called  Fermat  principle.  The  equations  of  the  eikonal 
and  bean  are  obtained  as  a consequence  of  this  principle  as  a result 
of  the  solution  to  variational  problem. 

In  the  case  of  parametric  media,  it  is  logical  to  require  the 
minimality  of  integral  (9),  that  determines  the  propagation  time  of 
the  recorded  phase  of  the  wave  between  two  points  of  space.  Then  it 
is  not  difficult  to  ascertain  that  the  equations  of  Euler  for  this 
variational  problem  in  accuracy  coincide  with  the  differential 
equation  of  beam.  In  other  words,  Fermat*  s principle  retains  its 


1 ’ 
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The  extremality  of  beans  nahes  it  possible  to  obtain  a 
convenient  fornula  for  an  instantaneous  signal  frequency  from  passage 
by  the  latter  of  path  s fron  the  starting  surface  on  which  we 
considered  it  nonochromatic  (with  frequency  wg) » i.e.,  « = ugfo  ♦ * 
on  Sq  (6  - the  certain  phase  shift).  By  definition  of  instantaneous 
frequency  at  the  aonent  of  tine  t 

"=(«)-= ""I®)/ 


I 


DOC  = 77218400 


PAGE 


Page  82. 


The  calculation  of  the  last  partial  derivative  (on  the  basis  of 
formula  (9))  substantially  is  facilitated  on  the  strength  of  equality 
to  *ero  of  first  variation  in  functional  (9)  with  the  attached  upper 
limit  and  the  "natural"  [4]  conditions  at  the  lower  end  of  the  path 
of  integration-  Finally  is  obtained  the  following  foraula: 

‘ = “«exp(j[S]rfsj.  (12> 


where  V*  - derivative  of  the  wave  propagation  velocity  in  terms  of 
tiie.  The  "departure"  of  frequency  in  the  parametric  medium  is  an 
integral  effect,  which  is  accumulated  along  the  entire  route  of  the 
propagation  of  beam.  But  with  the  limited  dimensions  of  parametric 
region  and  smallness  of  relative  changes  in  the  speed  (most  probable 
in  practical  sense  situation)  frequency  change  is  small. 
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§3.  Aaplitude  of  a lectronagnetic  field. 


Por  finding  the  amplitude  vectors  and  //o  we  have  the 


following  systea  of  equations,  which  ensues  fcoa  (1),  (5)  and  (6): 

It  dH«  II 

— * - ■■0* 

(13) 


rotE,+  f •- 


dt  e 


where  we  have  introduced  the  designation: 

. 1 dt  I d|t 

**£  — ’l  + TT  ■ dt'  “ in'  dt  • 


(14) 


The  first  of  values  (®£)  let  us  call  electrical,  and  the  second  (»«)  - 
peraeance  of  the  paraaetric  aediua. 

Otiliring  standard  transfer  nations  [4],  it  is  possible  to  pass 
froB  (13)  to  the  quadratic  relationship  which  characterizes  the 
energy  balance  in  terns  of  anplitude  vectors  (let  us  emphasize  the 
substantiality  of  the  latter)  : 

dl*  IE.  X H.I  (H.  • ^)  + I (e.-  f ) + ? (•,«:  + ».<)  = 0- 

(15) 

Fron  a aatheaatical  point  of  view  in  the  extracted  formula  «nter 
the  terns  only  of  first  (and  higher)  order  of  snallness  from  any  of 
three  low  paraaetars  of  our  problem. 
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Page  8 3. 


Let  as  coapose  the  scalar  functioa 

nuaerically  equal  to  the  energy  current  density  in  bean  direction. 
Taking  into  account  that  with  an  accuracy  to  terras  of  firs<--order 
saallness  the  vectors  and  satisfy  relationships  (0)  that 

deteraine  the  structure  of  the  field  of  plane  electromagnetic  wave, 
we  can  write  (15)  in  the  following  form  (for  the  wave,  which  goes  in 
direction  )‘. 

dlv(«l)  + l.  J + + = (15a) 

Equation  (15a)  is  distinguished  from  precise  formula  (15)  only 
by  terms  of  the  second  order  of  saallness  and  reflects  the  law  of 
conservation  of  energy,  which  takes  place  along  beam  tubes. 

Actually,  let  us  examine  the  dynamic  picture  of  beams,  which 
corresponds  to  the  recorded  starting  time  tor  and  let  us  isolate  beam 
tube  with  as  small  a cross  section  as  is  convenient  (Fig.  3).  Let  oj 
and  #2  of  i-ts  differentially  close  sections,  which  actually 

exist  in  times  t and  t ♦ dt  respectively  and  distant  from  each  other 
at  a distance  ds  = vdt.  Integrating  (15a)  at  moment  of  time  t in 
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connection  with  of  the  voluae  of  the  beaa  tube,  included  between 
these  sections,  utilizing  the  theoren  of  Gauss  and  retaining  only  the 
first-order  terns  of  snallness  with  respect  to  dt,  we  will  obtain 

lurfa  — j ud<i-^^^^(i->ds-\-<i^udzdsz=o^  (16) 


where 

the  total  flux  of  the  energy  through  a certain  section 
« of  bean  tube,  distant  at  a distance  s with  respect  to  the  bean  from 
starting  surface,  then  (16)  will  take  the  torn 

w(f.  s + rfs)  — w(^  + = 
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Fig.  3. 


Page  84. 

After  replacing  w(t,  s ♦ ds)  - v(t,  s)  by  we  will  oPtain 

dw-^-awds  = 0,  (16a) 

where  dw  is  the  tatal  differential  of  energy  flow  between  sections  oj 
and  taking  into  account  both  space  and  tiae  changes. 

Equation  (16a)  is  a differential  first-order  equation, 
describing  the  change  in  the  energy  flow  along  the  beam  tube.  It  is 
substantial,  that  it  is  derived  only  for  the  family  of  beams  united 
by  the  total  starting  time  to.  The  energy  flow  itself  is  giver,  by  the 
f oraula 

S 

- 

W(t,  s)  = «'o(/o,  0)e  “ , (17) 

where  w(to,  0)  - energy  flow  in  the  beginning  of  beam  tube,  and 
fading  is  determined  by  the  delaying  values  of  electrical 
conductivity  and  permeance  (14).  An  increase  in  the  permeability  of 
the  mediuB  leads  to  the  supplementary  energy  absorption  of  wave,  and 
a decrease  in  values  t and  p is  connected  with  the  release  of 
energy  by  the  medium. 
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For  tubes 
obtain  w = ud* 
beam  tube 


with  differentially  narrow  cross  section  do,  wp 
and  energy  current  density  in  certain  section  s 


u{t,  s)  = «„(/,.  0) 


rioCtr,  0) 
•to  (t.  S) 


(18) 


will 

ot  ’’ho 


In  formula  (18)  is  clearly  expressed  the  amplitude  reduction  or 
field  because  of  the  divergence  of  beams.  Since  it  is  obvious  tha^ 


u = 


(19) 


the  values  of  the  strength  of  fields  ^ and  no  easily  are  located 
through  (18). 


For  explaining  the  polarization  of  wave,  let  us  introduce  unit 
vectors  ^ and  /\  determining  at  each  point  of  space  the  sense  of 
the  vector  f,  and  Then  from  relationship  (6)  ensues  the 

connection  of  these  unit  vectors  with  the  unit  vector  of  beam  |* 

1 Xe  = -h:  lXh=e.  (20) 

Differentiating  the  extracted  relationship  with  respect  to  s,  we  will 
obtain 


V e 4- 1 y 


ds 


^ yh  + IX  -="* 


(21) 


1 


Assuming  the  polarization  of  the  wave  to  be  linear  (which  will 
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not  be  fundamental  limitation)  and  by  talcing  into  account  the  easily 
proven  fact  of  tha  retention  of  polarization  in  the  isotropic 
homogeneous  medium,  we  will  search  for  such  solutions  (21)  which 


satisfy  obvious  condition  ^ =0,if  ^ = 0- 


Page  85. 


Prom  a physical  point  of  view,  this  means  that  the  rotation  of 

vectors  ^ and  fi  along  beam  is  caused  only  by  bending  of  the 

latter,  but  formally  forces  from  all  possible  solutions  of  the  system 

Iv  — — — iv-**— 

‘ ^ as  ~ as  ’ as  ~ as 

to  select  only  the  zero  one.  Therefore  instead  of  the  complete 
equations  in  formula  (21)  we  can  examine  the  "shortened" 


relationsh  ips 


at  ah  at  ^ rfe 


Utilizing  (20),  we  pass  to  two  differential  equations: 

‘Oi_  ,/rfi  .]  at  .-at  \ 

^ (*•*’] . as  = ->  [as  • ••) ' (22) 

either  of  which  can  be  utilized  for  finding  the  polarization  of  wave 
with  known  radiation  picture. 


Specifically,  if  beam  lies  in  certain  plane  L,  then  from  (22)  it 
follows  that  the  electromagnetic  vector,  orthogonal  to  this  plane. 
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retains  its  direction.  In  this  case  there  is  no  need  to  integrate 
equations  (22),  and  it  suffices  to  separate  the  studied  field  into 
two  waves,  in  one  of  which  the  vector  is  orthogonal  to  L,  and  ir. 

the  other  ~ vector  • 

Equations  (22)  can  be  written  somewhat  otherwise;  talcing  into 
account  (11)  they  take  the  form 


^ = (22a) 


54.  Conditions  of  the  applicability  of  geometric-optical  observation. 


Thus,  we  obtained  the  following  picrture  of  the  propagation  of 
waves  in  a slightly  inhomogeneous  and  weakly- parametric  medium; 

1)  phase  distribution  of  wave  in  space  is  given  by  the 
generalized  equation  of  the  eikonal  (7) ; the  geometry  of  beams  is 
changed  in  the  course  of  time;  it  is  physically  reasonable  to  examine 
the  families  of  beams  which  correspond  to  one  starting  time; 

2)  the  energy,  which  flows  along  beam  tubes,  changes  its  value 
only  due  to  losses  in  the  medium  (the  latter  sometimes  can  be 
negative),  i.e.,  the  reverse  reflected  waves  do  not  appear. 


The  only  inaccuracy  in  our  reasonings  was  the  replacement  of 
formula  (15)  by  expression  (15a).  From  a physical  point  of  view,  this 
error  is  wholly  caused  by  the  difference  between  the  structure  of  our 
electromagnetic  field  and  the  structure  of  the  plane  of  wave  in  the 
homogeneous  medium. 

Page  86. 

For  plane  electromagnetic  waves,  as  is  well  known,  the  laws  of 
geometric  optics  are  precise. 

From  formal  point  of  view,  we  disregarded  terms  of  the  second 
order  of  smallness.  However,  we  have  not  yet  shewn  that  the 
satisfaction  of  inequalities  (2),  (3)  and  (4)  provides  the  accepted 

by  US  separation  of  values  by  order  of  smallness  (see  §1).  Now  this 
is  not  difficult  to  do,  since  we  have  formulas  (15)  and  (19),  that 
make  it  pciwible  to  produce  the  necessary  estimations  i. 


FOOTNOTE  *.  The  datailed  conducting  of  all  estimations  would  require 
much  space;  therefore  we  will  be  restricted  to  the  fundamental 
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1 

i 


(18)  and  (19).  ENDFOOTNOTE. 


For  example,  it  is  necessary  to  show  that 


OXi 


« 


IE„| 

X • 


(23) 


The  conditions  of  weak  heterogeneity,  weak  parametricity,  and 
smallness  of  losses  in  the  medium  assure  satisfaction  of  these 
inequalities,  but  with  one  supplementary  limitation,  namely: 

I V I ^ T ' it 

i.e.  we  must  examine  only  such  fields  for  which  a relative  change  ii. 
the  value  of  cross  section  d»  of  beam  tube  is  small  at  distances  of 
the  order  of  wavelength  and  at  times  on  the  order  of  the  period  of 
the  fluctuations  of  field.  A similar  requirement  will  be  carried  out, 
if  the  radii  of  curvature  of  phase  fronts  are  great  in  comparison 
with  wavelength.  In  other  words,  relationships  (24)  forbid  th®  use  of 

I‘*  geometric- optical  approximations  near  concentrated  sources  (necessary 

compulsorily  to  leave  into  their  wave  zone),  and  also  in  the  vicinity 

of  foci  and  caustic  surfaces  of  wave  field  (where  d®  ^ 0,  i.e., 

I the  beam  tubes  boundlessly  become  narrow).  Let  us  note  at  the  same 

^ time  that  geometric  optics  is  inapplicable  in  those  regions  of  space 

L (generally  speaking,  moving  with  time)  , 

I 

permeability  of  the  medium  (e  or  p) 

( 

I 


where  \ 4 «,  i.e.,  the 

vanishes. 


\ 


bL 


DOC  = 77218400 


PAGE 


The  second  ramarkable  fact  is  the  existence  of  the  integral 
inequalities,  which  ensure  evaluations  of  the  type  of  (23)  . For 
example,  the  condition  of  the  slowness  of  the  rotation  of  the  plane 
of  polarization,  which  is  reduced  on  the  strength  of  §3  to  the 
requirement  for  the  slow  rotation  of  beam,  taking  into  account  (11) 
takes  the  form 

0 ^ 

Page  87. 

Analogous  character  bears  the  inequality,  which  limits  the 
gradients  of  amplitude  between  adjacent  beam  tubes  and  requiring  that 
the  difference  in  fading  between  them  would  be  little: 

||£;I*'«T'  (»■» 

If  (26)  does  not  obtain,  then  it  is  not  possible  to  disregard  the 
effect  of  transverse  diffusion  [5]  between  beam  tubes. 

It  would  seei  that  inequalities  (25)  and  (26),  ensuring  the 
possibility  of  neglect  of  some  terms  in  the  differentia]  equations  of 
problem,  make  the  conditions  of  the  applicability  of  geometric  optics 
more  rigid  as  compared  with  the  case  of  the  nonconducting  and 
nonparametric  medium.  However,  this  is  not  so.  It  is  well  known  that 
the  error  in  solution,  caused  by  an  inaccuracy  in  the  equations,  can 


^ 1 
dt  <ii 


w 


DOC  = 77218400 


PAGE 


be  accunulated  during  considerable  changes  in  the  independent 
variables.  Actually,  if  we  instead  of  approxination  (15a)  writ<=‘  +he 
precise  relationship 

dlv(ul)  + l . ^+au  = /^(s,  ^),  (27) 

where  the  standing  in  right  side  function  F(s,  t)  - the  value  of  thf^ 
2nd  order  of  smallness,  is  egual  to 

F{s,  t)  = div  («1  — Eo  X Hnl  4“ 

-4-JLf^  ' -1  /7  ^^0  _ 1 1 / » , I 

I'  2 y i dt  2 f j 


to  evaluate  corrections  to  the  approximation  of  geometric  optics  on 
the  basis  of  formula  (27)  it  is  possible  to  write  the  integral 
relation  sh  ip 


«,(s,  t)do(s,  0 = J l/^(s')</»(s')lexpj-J’  l*l  r/s"j£(s'  + 

+ J (F(5')  do  (s'))  exp  j - ] 1»|  rfs"!  ds'.  (28) 


where  u is  the  difference  between  the  solutions  to  precise  (27)  and 
geometric- optical  (15a)  equations. 


To  ensure  the  relative  smallness  of  this  correction  (in  the 
general  case)  is  possible  only  under  the  condition  of  the  examination 
of  the  localized  heterogeneities  whose  linear  dimensions  D are 
commensurable  (but  not  greater!)  in  order  of  magnitude  with  the 
scales,  which  characterize  the  degree  of  heterogeneity  or 
parametricity  of  the  medium. 


PAGE 
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It  is  obvious  that  F(s,  t)  = 0(u  X/L^ ) , and  the  integral  in  formula 
(28)  can  be  estinated  as 


= 0^), 


where  D it  is  determined  by  the  length  of  that  part  of  the  beam  on 
which  is  essential  the  account  of  correction  F (s,  t) , i.e«,  D is  the 
dimension  of  heterogeneous  region.  The  smallness  of  ratio  XD/L 
provides  simultaneously  the  execution  of  integral  inequalities  (25) 
and  (26) , since  in  them  integrand  also  has  the  2nd  order  of 
smallness. 


Everything  said  in  this  section  makes  it  possible  to  assert  that 
the  fundamentally  new  limitations  for  applying  the  method  of 
geometric  optics  to  the  analysis  of  the  propagation  of 
electromagnetic  waves  in  the  parametric  medium  does  not  appear. 
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